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ABSTRACT 

We present high-dispersion spectroscopic data of the compact planetary nebula 
Vy 1-2, where high expansion velocities up to 100 km s -1 are found in the Ha, [N ii] 
and [O ill] emission lines. HST images reveal a bipolar structure. Vy 1-2 displays a 
bright ring-like structure with a size of 2.4" x3.2" and two faint bipolar lobes in the 
west-east direction. A faint pair of knots is also found, located almost symmetrically on 
opposite sides of the nebula at PA=305°. Furthermore, deep low-dispersion spectra are 
also presented and several emission lines are detected for the first time in this nebula, 
such as the doublet [Cl hi] 5517, 5537 A, [K iv] 6101 A, C ii 6461 A, the doublet C iv 
5801, 5812 A. By comparison with the solar abundances, we find enhanced N, depleted 
C and solar O. The central star must have experienced the hot bottom burning (CN- 
cycle) during the 2 nd dredge-up phase, implying a progenitor star of M^3 M@. The 
very low C/O and N/O abundance ratios suggest a likely post-common envelope close 
binary system. 

A simple spherically symmetric geometry with either a blackbody or a H-deficient 
stellar atmosphere model is not able to reproduce the ionisation structure of Vy 1-2. 
The effective temperature and luminosity of its central star indicate a young nebula 
located at a distance of ~9.7 kpc with an age of ^3500 years. The detection of stellar 
emission lines, C II 6461 A, the doublet C ivAA 5801, 5812 and O III 5592 A, emitted 
from a H-deficient star, indicates the presence of a late-type Wolf Rayet or a WEL 
type central star. 

Key words: ISM: jets and outflows - ISM: abundances - planetary nebula: individ¬ 
ual: Vy 1-2 - stars: Wolf-Rayet 


1 INTRODUCTION 

Planetary nebulae (PNe) are the consequent result of the in¬ 
teraction of stellar winds ejected from low-to-intermediate 
mass stars (1-8 M@). A slow and dense wind (10 -4 Mg), 
emerges from an Asymptotic Giant Branch (AGB) star, is 
swept up into a dense spherically symmetric shell by a sec¬ 
ond, faster and tenuous wind expelled during the post-AGB 
phase (interacting stellar wind model; Kwok et al. 1978). 
This model can reasonably explain the formation and evo¬ 
lution of spherical PNe. However, most post-AGB stars and 
PNe show highly asymmetric circumstellar envelopes and 
complex morphologies (Manchado et al. 2006, Meixner et 
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al. 1999, Parker et al. 2006, Boumis et al. 2003, 2006, Mis- 
zalski et al. 2008, Lagadec et al. 2011, Sahai 2011, Sabin et 
al. 2014). 

The deviation of PNe from spherically symmetric struc¬ 
tures has its origin in an enhanced equatorial-to-polar den¬ 
sity gradient in the AGB circumstellar envelopes resulting 
in the formation of aspherical and more complex nebulae. 
This equatorial density enhancement is interpreted as the 
result of common envelope evolution of close binary systems 
(Soker & Livio 1994, Nordhaus & Blackman 2006, De Marco 
2009). Although, most of the post-AGB and PNe show com¬ 
plex morphologies, only a small number of them are known 
to contain a binary system (Miszalski et al. 2009, De Marco 
2009). Recent work by Jones et al. (2014) confirmed the 
binary nature of Hen 2-11’s nucleus, providing further evi- 
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DEC 



24.0 23.5 17:54:23.0 22.5 22.0 


R.A. 

Figure 1. Slit positions of the low—(1) and high-(2, 3, 4 and 5) dispersion spectra superimposed upon the high-resolution HST WFPC3 
narrow-band [O III] image (ID: 11657; PI: Stanghellini). Left panel reveals the presence of a pair of knots at P.A.=305°. Right panel is 
in logarithmic scale and illustrates the ellipsoidal structure of the nebula with extended filamentary structures and the presence of two 
fainter lobes at P.A. ~90°. 


dence of a positive link between binary nuclei and aspherical 
PNe. 

The compact planetary nebula Vy 1-2 (PNG 
053.3+24.0), discovered by Vyssotsky (1942), had first been 
described as an elliptical PN (“The I AC Morphological Cat¬ 
alogue of Northern Galactic Planetary Nebulae”; Manchado 
et al. 1996), but recent high angular resolution HST images 
revealed a more complex structure (Sahai et al. 2011, Shaw 
et al. in preparation). According to Winberger’s catalogue of 
expansion velocities of Galactic PNe, Vy 1 -2 shows V[o in] 
> Vjn ii] (Weinberger 1989). This result is inconsistent with 
a simple homologous expansion law. A similar velocity trend 
has also been found for few other PNe, (e.g. BD +30°3639, 
Akras & Steffen 2012, Bryce & Mellema 1999; see also Med¬ 
ina et al. 2006 for more objects). A detailed study of the 
velocity field as a function of distance from the central stars 
in PNe has revealed a more complex expansion law than the 
simple homologous-law, probably associated with the fast 
wind from the central stars (Gesicki & Zijlstra 2003). In¬ 
terestingly, the central stars of this rare group of PNe have 
been identified as Wolf Rayet type ([WR]) or weak emission 
line stars (WELS) (Medina et al. 2006). Unfortunately, the 
spectral type of the central star of Vy 1-2 is still unknown. 
Barker (1978a,b), and more recently Wesson et al. (2005), as 
well as Stanghellini et al. (2006) studied the ionisation struc¬ 
ture and chemical composition of this nebula. A comparison 
between these studies shows a significant discrepancy in the 
ionic and total chemical abundances (see Tables 5 and 6). 
This may be related to the different values of T e and N e 
used by these authors. In particular, Wesson et al. (2005) 
estimated N e and T e using diagnostics line ratios from low 
to high excitation lines (see Table 4). The high difference in 
N e between [S n] and [O n] cannot be explained since they 
have similar ionisation potentials. Barker (1978a) gave an 
average value of 5000 cm -3 using different diagnostics from 
low excitation lines ([S n] and [O n]). Therefore, the low N e 


in [S n] of 1160 cm 3 derived by Wesson et al. (2005), may 
not be correct. 

The logarithmic H/i flux (erg s _1 cnh 2 ) was estimated 
to be -11.52+0.01 (Barker 1978a; hereafter Bar78a) and - 
11.53 (Wesson et al. 2005; hereafter WLB05). Vy 1-2 is also 
a radio source with 1.4 GHz flux equal to 12.4 + 0.9 mjy and 
extinction coefficient c(radio)=0.1, which is slightly lower 
than the optical value from WLB05 (c(H/3)=0.139). The co¬ 
ordinates of Vy 1-2 are given at RA 17 h 54 m 23 s and Dec. 
27° 59' 58" for the optical, and RA 17 h 54 m 23 a .49 and Dec. 
+27° 59 ; 58".5 for the radio position (Condon & Kaplan 
1998). 

The first estimation of the effective temperature of the 
central star of Vy 1-2 was by Gurzadyan (1988), 124 kK and 
66 kK using the He II / H f3 and [O ill]/ He II line ratios, re¬ 
spectively. More recently, Stangellini et al. (2002) calculated 
T eff = 119+25 kK and log(L/L 0 ) = 3.275+0.285, for a dis¬ 
tance of 7.24 ± 0.13 kpc, whereas Phillips (2003) estimated 
the He i and He ii Zanstra temperatures of 75.4 ± 10.2 kK 
and 99.8 ± 4.3 kK, respectively. Its distance is uncertain. It 
has been estimated to range from 4.38 kpc to 7.6 kpc: 4.7 kpc 
(Maciel 1984); 5.7 kpc (Sabbadin et al. 1984); 4.38 kpc 
(Phillips 1998) and 7.6 kpc (Cahn & Kaler 1971, Calm et al. 
1992). However, all the aforementioned authors used a radius 
of 2.4" for their distance estimation obtained from Cahn and 
Kaler (1971). No high-angular resolution of Vy 1-2 had been 
obtained until 1996 (Manchado et al. 1996). More recently, 
Shaw et al. (in preparation) obtained the high-angular HST 
images that we present here. The value of 2.4" is almost ~1.6 
times larger than the radius find in this work using the HST 
images (see below). This implies that the previous reported 
distances of Vy 1-2 are likely underestimated. 

The large differences in chemical abundances and ef¬ 
fective temperature inspired us to study this nebula and 
provide more accurate determinations of stellar and nebu¬ 
lar properties. We present high-resolution HST images of 
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Table 1 . Observation log 


Slit 

Dispersion 

Date 

Filter 


wavelength (A) 

FWHM (A) 

P.A. (°) 

Exp. times (sec) 

1 

low 

2011 Sep 04 

- 


1.3m Skinakas Telescope 
4700-6800 

- 

0 

3600 and 300x3 

2 

high 

2011 May 15 

Ho+[N 

ii] 

2.1m SPM Telescope + MES 
6550 

90 

90 

1800 

2 

high 

2011 May 15 

[O m] 


5007 

60 

90 

1800 

3 

high 

2012 Mar 30 

Ha+[N 

ii] 

6550 

90 

0 

1200 

4 

high 

2012 Apr 01 

Ha+[N 

ii] 

6550 

90 

55 

1200 

5 

high 

2012 Apr 01 

Ha+[N 

ii] 

6550 

90 

310 

1800 

5 

high 

2012 Apr 01 

[O in] 


5007 

60 

310 

1800 

6“ 

high 

2002 Jul 23 

Ha+[N 

ii] 

6550 

90 

310 

1800 


imaging 

2010 Feb 08 

F502N 


HST+WFPC3 

5010.0 

65.0 


600 

- 

imaging 

2010 Feb 08 

F200LP 


4883.0 

5702.2 

- 

11 

- 

imaging 

2010 Feb 08 

F350LP 


5846.0 

4758.0 

- 

11 

- 

imaging 

2010 Feb 08 

F814W 


8024.0 

1536.0 

- 

23 


a It was drawn from the SPM Kinematic Catalogue of Planetary Nebulae 


Vy 1-2 in the [O ill] emission line and broad-band filter 
F200LP, as well as new high- and low-dispersion spectra 
in order to provide further insight into the kinematic struc¬ 
ture and chemical composition of this compact nebula. In 
Section 2, we describe the observations and data reduction 
procedures. In Section 3, we present the results from the 
morpho-kinematic and chemical analysis. We also present 
the physical properties and chemical abundances derived 
from the ionisation correction factor (ICF; Kingsburgh & 
Barlow 1994) and from a photo-ionisation model performed 
with Cloudy (Ferland et al. 2013). Our results are discussed 
in Section 4 and conclusions are summed up in Section 5. 


2 OBSERVATIONS 
2.1 HST imaging 

High-resolution HST images of Vy 1-2 were drawn from the 
Mikulski archive for space telescopes (Shaw et al. in prepa¬ 
ration). They were taken with the Wide Field Planetary 
Camera 3 (WFPC3) on February 8 th 2010 (GO program 
11657, PI: Stanghellini L.). The narrow-band F502N filter 
(A=5009.8 A) was used to isolate the [O hi] 5007 A emission 
line and determine the morphology of Vy 1 2 (Fig. 1), while 
the broad filters F200LP, F350LP and F814W (“I”, see Ta¬ 
ble 1) were obtained in order to probe the presence of a cold 
binary companion. The exposure times were 600, 11, 11, and 
23 sec, respectively. A faint outer shell or halo as well as a 
pair of knots can be discerned in the [O ill] image (Fig. 1, 
left panel), while the F200LP (Fig 1, right panel) image is 
used to highlight the inner nebular components: the bright 
elliptical ring-like structure and the faint lobes. 

We carefully investigate whether the outer shell/halo 
is a real nebular component or the result of scattered light 
from the bright inner nebula in the HST optics. In Figure 2, 


left panel, we display the radial profile of a theoretical PSF0 
(green solid line), of two held stars (observed with the same 
camera and filter; red and blue solid lines) and of Vy 1- 
2 (yellow dashed line). The radial profiles between the two 
held stars and the theoretical PSF agree reasonably well. A 
discrepancy with the modelled PSFs can be seen at around 
1.5-2”, though it is negligible compared to the Vy 1-2’s pro¬ 
file. On the other hand, Vy 1-2 shows a maximum brightness 
at r=1.2" (at the position of the ring-like structure), while 
for r>2.2-2.3” (indicated by a vertical arrow in Fig. 2) shows 
a similar steep decline for the stars and theoretical PSF but 
two magnitude brighter. The theoretical [O ill] PSF was 
used to convolve the HST image using Leon Lucy’s algo¬ 
rithm in STSDAS. The azimuthally-averaged prohle of the 
convolved image is presented in Figure 2, right panel. The 
bump at r=2.2-2.3” is indicated by a vertical arrow. For 
r>2.3”, the prohle is a constant function, which means that 
the outer shell/halo is seen due to the PSF scattering. In par¬ 
ticular, for a power-law density distribution of p(r)ocr _1 , a 
r -2 surface brightness is expected (Sahai et al. 1999), which 
we do not see here. Nevertheless, some contribution of scat¬ 
tering light from dust can not be excluded, given that Vy 1-2 
has been found to be a strong IR emitter. 

2.2 High dispersion longslit spectroscopy 

High-dispersion, long-slit spectroscopic data of Vy 1 -2 in 
Ha+[N n] and [O ill] were obtained using the Manchester 
echelle spectrometer, (MES-SPM; Meaburn et al. 2003) on 
the 2.1 m telescope at the Observatorio Astronomico of San 
Pedro Martir in Baja California, Mexico, in its f/7.5 con¬ 
figuration. The observing run took place on May 2011 and 
April 2012 (see Table 1). The MES-SPM was equipped with 
a Marconi CCD detector with 2048x2048 square pixels and 
each side of 13.5 /im (=0.176” pixel -1 ). A 90 A bandwidth 

1 The program Tiny Tim developed by Krist et al. (2011), was 
used to create the theoretical PSF for the [O III] image. 
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Figure 2. Left panel: surface brightness profile of [O III] emission in Vy 1-2 averaged over all angles (yellow dashed line). The bump at 
r~2.2 "—2.3 ", indicated by an arrow, corresponds to the contribution of the faint lobes along the east—west direction. The red and blue 
solid lines represent the radial brightness of two stars observed with the same camera and filter while the green solid line corresponds to 
the profile of a theoretical HST point spread function (PSF) created using the program Tiny Tim. Right panel: surface brightness profile 
of the convolved [O III] image averaged over all angles, using the theoretical PSF and Leon Lucy’s algorithm. 


filter was used to isolate the 87 th order containing the Ha 
and [N n] AA 6548, 6584 A, and the 113 t/l order contain¬ 
ing the [O ill] 5007 A emission lines. Two-by-two binning 
was employed in both the spatial and spectral directions (a 
three-by-three binning was employed in both directions for 
the observing run in May). Consequently, 1024 increments, 
each 0 / . , 352 long, gave a projected slit length of 5(12 on the 
sky. The slit used was 150 fim wide (= 11 km s -1 and l"9). 
During the observations the slit was oriented at different po¬ 
sition angles. The integration times were 1800 sec, 1500 sec 
and 1200 sec (see Table 1 for details). In order to avoid any 
contamination from the bright central star, the slits were 
slightly offset from it by ~0.5 ", covering only the bright 
ring-like structure. 

A faint continuum could be discerned within the very 
short spectral range of our single-order echelle spectra. The 
wavelength calibration was performed using a Th/Ar cali¬ 
bration lamp to an accuracy of ±1 km s -1 when converted 
to radial velocity. The data reduction was performed using 
the standard IRAF 0 and STARLINK routines. Individ¬ 
ual images were bias subtracted and cosmic rays removed. 
All spectra were also calibrated to heliocentric velocity. The 
position-velocity (PV) diagrams and line profiles for the Ha, 
[N n] and [O III], and for all slit positions are presented in 
Figures 4-8 and Figures 9-12, respectively. 



Figure 3. Schematic diagrams of Vy 1—2 using the code SHAPE. 
Panels a) to c) show the mesh SHAPE images for three differ¬ 
ent orientations right, front and top views, respectively. Panel d) 
illustrates the final rendered 2-D image. 


2.3 Low dispersion spectroscopy 

Low-dispersion spectra of Vy 1-2 were obtained with the 

1.3 m telescope at Skinakas Observatory in Crete, Greece, 
on September 4 th , 2011. A 1300 line mm -1 grating was used 

2 IRAF (Image Reduction and Analysis Facility) is distributed 
by the National Optical Astronomy Observatory, which is oper¬ 
ated by he Association of Universities for Research in Astronomy 
(AURA) Inc., under cooperative agreement with the National Sci¬ 
ence Foundation. 


in conjunction with a 2000 x 800 SITe CCD (15 ^tm pixel), 
resulting in a scale of 1 A pixel -1 and covering the wave¬ 
length range of 4700-6800 A. The slit width used was 7.7" 
with a slit length of 7.9 ' and oriented in the north-south di¬ 
rection (Slit position 1). Two different exposure times were 
taken in order to identify both the bright and faint emis¬ 
sion lines: 3600 sec (saturated in Ha and [O ill] 5007 A 
lines, but detecting the much fainter lines) and 3 x 300 sec. 
The spectro-photometric standard stars HR5501, HR7596, 
HR9087, HR718 and HR7950 (Hamuy et al. 1992) were ob- 
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Figure 4. Hq (upper panel), [N ii] 6584 A (middle panel) and [O ill] 5007 A (lower panel) position-velocity (PV) diagrams of Vy 1—2 
shown at two different contrast levels the for slit position 2 (P.A.=90°). 
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Figure 5. Ha (upper panel) and [N n] 6584 A (lower panel) position-velocity (PV) diagrams of Vy 1-2 shown at two different contrast 
levels for the slit position 3 (P.A.=0°). 
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Figure 6. Ha (upper panel) and [N n] 6584 A (lower panel) position-velocity (PV) diagrams of Vy 1-2 shown at two different contrast 
levels the for slit position 4 (P.A.=55°). 
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Figure 7. Ha (upper panel), [N ii] 6584 A (middle panel) and [O ill] 5007 (lower panel) position-velocity (PV) diagrams of Vy 1-2 
shown at two different contrast levels for the slit position 5 (P.A.=310°). 


© 2013 RAS, MNRAS 000, [I]-?? 



Evidence for a Wolf-Rayet or WEL-type binary nucleus in the bipolar Planetary Nebula Vy 1-2 7 



0 * 

-40(1 -300 -200 -100 0 100 200 -400 -300 -200 -100 0 100 200 


Vhel i kill s ' 1 ) V H el (kins' 1 ) 

Figure 8. Ha (upper panel) and [N II] 6584 A (lower panel) position-velocity (PV) diagrams of Vy 1-2 shown at two different contrast 
levels for the slit position 6 (P.A.=310°), drawn from the SPM catalogue (Lopez et al. 2012). 


served in order to flux calibrate the spectra, while the fluxes 
were also corrected for atmospheric extinction and interstel¬ 
lar reddening. 

The data reduction was performed using standard IRAF 
routines. Individual images were bias subtracted and flat- 
field corrected using a series of twilight flat-frames. The sky- 
background was also subtracted from each individual frame 
while no dark frames were used as the dark current was 
negligible. Finally, the wavelength calibration was performed 
using a Fe-He-Ne-Ar calibration lamp. 

3 RESULTS 
3.1 Morphology 

Figure 1 displays the high-resolution HST [O hi] (left panel) 
and ultraviolet (F200LP, right panel) images of Vy 1-2, at 
two different scales, in order to highlight the fainter nebular 
components. The left panel displays a faint outer shell/halo 
and two much fainter knots at a position angle (hereafter 
P.A.) of 305°, located at 11.5” (knot S) and 14.5” (knot N) 
from the central star. The right panel shows an enlargement 
of the nebula, revealing a bright elliptical ring-like structure 
with a minor and major axis of 2.4” and 3.2”, respectively. 
This implies a mean radius of ~1.4”. This value is ~1.6 
times lower than the value used previously for the distance 
estimation obtained from ground-based observations. Thus, 
we recalculated the distance of Vy 1-2 using our more ac¬ 
curate value of radius, and we obtain a distance between 
7.0 kpc and 14.5 kpc. These distances clearly indicate that 
Vy 1-2 is a distant nebula. An average distance of 9.7 kpc 
is adopted here. 

The ring-like structure has a non-uniform emission 
along its perimeter, with an apparent filamentary structure, 
which indicates a strong interaction of stellar winds in ad¬ 
dition of an inner cavity. Its major axis is oriented towards 
to a P.A.=320°, 15° rotated from the knots’ direction. This 


may be associated with a possible precession of the rota¬ 
tion axis of a single central star, or the orbital plane of a 
binary star at the centre of Vy 1-2. Two faint lobes are 
also apparent in the east-west direction (P.A. ~90°). De¬ 
spite their very low surface brightness compared to the cen¬ 
tral ring-like structure, their morphology suggests a bipolar 
structure, like the Homunculus nebula and Hb 5 PN. The 
eastern and western bipolar lobes are found to be blue- and 
red-shifted, respectively, while their different sizes imply a 
tilted nebula towards the line of sight. In Figure 3, we illus¬ 
trate a schematic 3-D diagram of Vy 1-2 for different points 
of view, using the code SHAPE (Steffen et al. 2011). This 
simple model reproduces fairly well the main components 
of Vy 1-2, such as the ring—like structure and the bipolar 
lobes as well as the orientation of the nebula (inclination and 
position angles). The inclination angle of the model is 13°, 
while the ring-like structure is best described with an ellipse. 
Any attempt to reproduce the ring-like structure assuming 
a round shape, resulted in a very high inclination angle of 
~50° that we do not get from the PV diagrams. Our 3-D 
structure of Vy 12 is only one of many other solutions that 
can be obtained just by varying the size of bipolar lobes and 
the orientation accordingly. 

Moreover, in order to investigate the possible pres¬ 
ence of a cold binary companion, the F200LP, F350LP and 
F814W broad-band images were studied. With the available 
data, we found no evidence of a nearby companion. If the 
central star is a binary system, then the projected separation 
must be smaller than 0.08”. 


3.2 Kinematics 

In Figures 4 to 8 we present the observed Ha 6563 A, [N ii] 
6584 A and [O ill] 5007 A position-velocity (PV) diagrams 
along the slit positions 2 to 6. All slit positions are shown 
overlaid on the HST [O ill] image (Fig. 1, left panel). Besides 
the common bright Ha, [N ii] and [O ill] emission lines, the 
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Figure 9. The Ha, [N II] 6584 A and [O III] 5007 A line profiles for slit position 2 (P.A.=90°) are shown in the left, middle and right 
panels, respectively. The blue and green lines correspond to the observed line profile and a single Gaussian best fit, respectively. 




Figure 10. The Ha and [N n] 6584 A line profiles for slit position 3 (P.A.=0°) are shown in the left and right panels, respectively. The 
blue and green lines correspond to the observed line profile and the single Gaussian best fit, respectively. 




Figure 11. The Ha and [N n] 6584 A line profiles for slit position 4 (P.A.=55°) are shown in the left and right panels, respectively. The 
blue and green lines correspond to the observed line profile and the single Gaussian best fit, respectively. 





Figure 12. The Ha, [N ii] 6584 A and [O ill] 5007 A line profiles for slit position 5 (P.A.=310°) are shown in the left, middle and right 
panels, respectively. The blue and green lines correspond to the observed line profile and the single Gaussian best fit, respectively. 
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Table 2. Expansion velocities of Vy 1-2 


slit 

RA. 

(°) 

V(H a) 
km s —1 

V([N ii]) 
km s 1 

V([0 m]) 
km s 1 

V(He n) 
km s —1 

low velocity component 



2 

90 

19 

9 

18 

25 

3 

0 

18 

8 

- 

- 

4 

55 

19 

9 

- 

25 

5 

310 

19 

9 

18 

23 

6 

310 

19 

8 

- 

23 

high velocity component 



2 

90 

90 

40 

90 

- 

3 

0 

100 

40 

- 

- 

4 

55 

100 

40 

- 

- 

5 

310 

95 

90 

90 

- 

6 

310 

100 

90 

- 

- 


He ii 6560 A and He I 5015 A lines, as well as the much 
weaker lines, C n 6581.2 A, C n 6578 A, [N n] 4996.8 A, 
[Fe i] + [Fe Ii] 5001.9 A or [N n] 5002.7 A, [N n] 5005.2 A 
and [Fe ii] 5018.2 or He I 5017.2 A are also present (see Figs 
4 to 8). 

By conducting a literature search for these emission 
lines in well known objects, we found a small number of 
PNe (e.g NGC 6543, Hyung et al. 2000; M 2-9, Lee H.-W. 
et al. 2001) and symbiotic stars (RR Telescopii and He 2- 
106, Lee H.-W. et al. 2001) that exhibit the 6581 A emission 
line. However, those objects also exhibit the 6545 A line, 
which usually is attributed to the Raman scattered He II 
line that is not detected in our data. Both the 6578 A and 
6581 A lines have also been detected in the symbiotic star 
AE Circinus, and identified as carbon lines (Mennickent et 
al. 2008). The 5002 A line has been detected in the J320 
(PNG190.317.7) nebula but no identification was provided 
(Harman et al. 2004). 

The PV diagrams clearly show the presence of a high¬ 
speed component (wide wings) with velocities up to 90- 
100 km s _1 , and a much slower component ~20 km s _1 
(Gaussian fit). The Ha 6563 A, [N n] 6584 A and [O hi] 
5007 A emission line profiles, in Fig. 9 to 12, are from the slit 
positions 2 to 5, respectively. Single Gaussian fit to these line 
profiles gives turbulent widths (FWHM) of 19±3 km s' 1 
in the Ha, between (8-9)±2 km s' 1 in the [N ii] , and 
18±3 km s' 1 in the [O hi] emission (see Table 2). These 
expansion velocities correspond to the bright, inner ring like 
structure. Two separated blobs can be discerned in the [N ii] 
6584 A PV diagram from the slit position 4 (Fig. 6; right, 
lower panel), similar to those of a pole-on toroidal or ring¬ 
like structure, whilst high expansion velocities are not ap¬ 
parent, since this slit position covers only a small part of the 
bipolar lobes. All the widths are corrected: for the 20 km s -1 
instrumental width; the 21.4 km s -1 thermal width of the 
Ha line at 10000 K (5.8 km s' x and 5.4 km s _1 of the [N Ii] 
and [O ill] lines respectively); and the fine structural broad¬ 
ening of 6.4 km s' 1 of the Ha line. Note that, unlike the 
Ha line, the [N Ii] and [O ill] lines have no fine structural 
components. The He II expansion velocity is estimated (23- 
25)±4 km s -1 from the He II 6560 A recombination line, and 
is corrected for the instrumental broadening of 20 km s' 1 , 


the thermal broadening of 10.7 km s' 1 (10000 K), and 
the fine structural broadening from this line at ~6 km s' 1 
(Meaburn et al. 2005). 

Additional to the previous low velocity components, the 
line profiles of Vy 1-2 also present broad wings (Figs. 9-12). 
This continuous high-velocity component likely corresponds 
to the expansion of the faint bipolar lobes seen in the HST 
images. V[o nil is between 90±4 km s _1 and 100±4 km s' 1 , 
measured at the base of line profiles for all the P.A. The 
[N ii] emission line shows such high expansion velocity of 
90-100 km s' 1 only along the P.A.=310°, whilst for the rest 
of the P.A.s, V[n ii] is significantly lower than V[o mj. This 
is in agreement with the Winberger’s catalogue (Weinberger 
1989). A similar velocity trend, V[o in] > V[ N n], has also 
been reported for a small group of PNe (see §4.1). Regarding 
the HST images and PV diagrams, we conclude that Vy 1- 
2 is seen almost pole-on with an inclination angle of 10° 
±2° to the line of sight. The average systemic velocity of 
the nebula is 14 ys =-85 km s' 1 , with a standard deviation 
of 4 km s' 1 , in good agreement with the previous value of 
-82.4±3.9 km s' 1 quoted by Schneider et al. (1983). 

Figure 8 shows the PV diagram obtained along 
P.A.=305° in Ha+[N ii], drawn from the recently released 
SPM Kinematic Catalogue of Galactic Planetary Nebu¬ 
lae (Lopez et al. 2012). This spectrum covers, apart from 
the bright ring-like structure, the faint knot N (Fig. 1). 
The Ha line profile of this knot is centred on Vhel=- 
80±4 km s' 1 . This central value should be compared to 
the systemic heliocentric radial velocity of the whole neb¬ 
ula, V ays =-85±4 km s' 1 . 


3.3 Ionisation structure 

3.3.1 Nebular parameters and ICF abundances 

The panels in Figure 13 display the low-dispersion spectra of 
Vy 1-2 that cover the wavelength range of 4700-6800 A for 
exposure times of 3600 (upper panel) and 300 sec (lower pan¬ 
els), respectively. Several emission lines from low-to-high 
and very high ionisation states are present and some are de¬ 
tected for the first time in this nebula. Besides the common 
nebular lines, such as Ha, H/3, He i, He II, [O ill], [N ii] and 
[S ii] lines, the [Cl hi] 5517 and 5537 A, C ii 6461 A, C iv 
5801 and 5812 A, [Si n] 5040 A, [K iv] 6101 A and [N Ii] 
5679 A emission lines are also detected. 

It is worth mentioning here that the emission line at 
5679 A has been identified either as the [N n] line by Liu et 
al. (2000; NGC 6135, 2001; Ml-42 and M2-36) and Ercolano 
et al. (2004; NGC 1501), or as [Fe vi] line by Hyung & Aller 
(1997). However, the latter classification requires, a very hot 
star of T e ff ^140 kK, which is not in agreement with the 
values previously published (T e g=119 kK, Stangellini et al. 
2002; 75.4 kK< T eff <99 kK, Phillips 2003). 

All very faint emission lines presented here are only de¬ 
tected in our deep spectrum (t eX p=3600 sec). Notice that 
WLB05 did not detect any of these lines due to their shal¬ 
low spectroscopy. 

In Table 3, we list the optical nebular emission line 
fluxes, corrected for atmospheric extinction and interstellar 
reddening. All line fluxes are normalized to F(H/3)=100. The 
interstellar extinction c(H/3) was derived from the Balmer 
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Table 3. Observed and modelled emission line intensities for 
Vy 1-2, corrected for interstellar extinction. The lines outside 
the wavelength range 4700-6800 A were drawn from WLB05. 


Line(A) 

1(A) 

error 

model 

WLB05 

[O n] 3726 

— 

— 

17.69 

20.04 

[O n] 3729 

- 

- 

9.57 

9.97 

[Ne hi] 3869 

- 

- 

98.86 

94.55 

[Ne hi] 3968 

- 

- 

29.81 

29.91 

He I 4026 

- 

- 

2.05 

2.19 

[S n] 4070 

- 

- 

1.49 

1.21 

Hq 4100 

- 

- 

25.89 

24.44 

H<5 4340 

- 

- 

46.81 

45.30 

[O hi] 4363 

- 

- 

9.12 

8.51 

He i 4471 

- 

- 

4.39 

4.40 

He ii 4686 

- 

- 

26.85 

31.00 

He I 4712 

- 

- 

0.54 

0.52 

[Ar iv] 4712 

3.89 

0.43 

3.38 

4.29 

[Ar IV] 4740 

3.97 

0.48 

3.72 

4.22 

H/3 4861 

100.00 

4 

100.00 

100.0 

He I 4924 

1.38 

0.34 

1.18 

1.15 

[O hi] 4959 

434.85 

13 

497.21 

443.10 

[N ii] 4997*’ 

- 

- 

- 

- 

[Fe ii] 5002 6 

- 

- 

- 

- 

[N ii] 5003 6 

- 

- 

- 

- 

[N ii] SOOS 6 

- 

- 

- 

- 

[O hi] 5007 

1287.06 

38 

1496.59 

1337.0 

He i 5015 b 

- 

- 

- 

- 

He i 5017 6 

- 

— 

— 

— 

[Fe ii] 5018 6 

- 

- 

- 

- 

[Si ii] 5040“ 

0.25 

0.09 

- 

- 

[N i] 5200 

0.22 

0.07 

0.00 

- 

He II 5411 

2.15 

0.19 

1.95 

- 

[Cl hi] 5517 

0.56 

0.11 

0.52 

- 

[Cl hi] 5537 

0.64 

0.16 

0.61 

- 

[O i] 5577“ 

0.05 

0.03 

0.001 

- 

[N ii] 5678 

0.13 

0.04 

0.00 

- 

[N ii] 5755 

0.53 

0.11 

0.42 

- 

He i 5876 

13.16 

0.71 

13.11 

13.10 

[K iv] 6101“ 

0.23 

0.06 

- 

- 

He ii 6119 

0.11 

0.04 

0.06 

- 

He ii 6172 

0.08 

0.04 

0.06 

- 

He II 6234 

0.08 

0.04 

0.07 

- 

[O i] 6300 

2.02 

0.14 

0.06 

1.70 

[S m] 6312 

1.45 

0.15 

5.09 

1.30 

[O i] 6363 

0.69 

0.09 

0.02 

0.71 

He ii 6405 

0.11 

0.04 

0.11 

- 

[Ar v] 6435 

0.17 

0.06 

0.43 

- 

He ii 6525 

0.17 

0.06 

0.13 

0.5 

[N ii] 6548 

8.42 

0.45 

8.75 

9.16 

He ii 6560 b 

- 

- 

- 

- 

Ho 6563 

284.62 

10 

287.45 

284.1 

C ii 6578 6 

- 

- 

- 

- 

C ii 6581 6 

- 

- 

- 

- 

[N ii] 6584 

26.23 

1.32 

25.81 

22.68 

He i 6678 

3.79 

0.31 

3.59 

3.89 

[S ii] 6716 

2.49 

0.28 

2.61 

2.04 

[S ii] 6731 

4.13 

0.32 

4.69 

2.47 

[Ar v] 7005 

- 

- 

0.95 

0.61 

He i 7065 

- 

- 

4.39 

3.96 

[Ar in] 7135 

- 

- 

12.51 

13.83 

He i 7281 

— 

— 

0.70 

0.50 

[Ar hi] 7751 

- 

- 

3.01 

2.61 

logF(H/}) 

-11.56 

— 

-11.56 

-11.53 


Table 3. Continued 


Line(A) 

1(A) 

error 

model 

WLB05 

[O III] 5007/4363 

- 

— 

164 

157 

[N ii] 6581/5755 

49.5 

10.4 

61.4 

- 

[O i] 6300+6363/5577 

54.2 

28.1 

80.0 

- 

[O ii] 3727/3729 

- 

- 

1.85 

2.01 

[S ii] 6716/6731 

0.60 

0.12 

0.56 

0.83 

[Ar iv] 4712/4740 

0.98 

0.27 

0.92 

1.01 

[Cl hi] 5517/5537 

0.88 

0.36 

0.85 

- 

J 1.25“ 

3.5 d 

0.2 d 

2.4 

— 

H 1.65 

2.0 

0.1 

2.3 

- 

K 2.15 

2.6 

0.1 

2.8 

- 

WISE 3.35 

3.4 

0.1 

2.5 

- 

WISE 4.60 

5.6 

0.1 

2.2 

- 

AKARI 9.00 

104 

21 

65 

- 

WISE 11.6 

164. 

21 

110 

- 

IRAS 12.0 

113 

22 

211 

- 

AKARI 18.0 

786 

21 

893 

- 

WISE 22.1 

1319 

15 

1804 

- 

IRAS 25.0 

1647 

83 

2361 

- 

IRAS 60.0 

2078 

104 

2122 

- 

AKARI 65.0 

1704 

_e 

1795 

- 

AKARI 90.0 

1184 

58 

796 

— 

IRAS 100.0 

1094 

219 

592 

- 

AKARI 140.0 

697 

_e 

218 

- 


a These lines were not used to constrain the model. 
b Emission lines detected in the high-dispersion spectra 
c In /mi, d In mjy 

e No errors are provided for these measurements since the flux 
quality indicator is 1. 


Ha/H/3 ratio (Eq. 1), using the interstellar extinction law 
by Fitzpatrick (1999) and R v = 3.1, 


<H0) 


0.348 


log 


F(Ha)/F(H/3) 

2~85 


(1) 


where the 0.348 is the relative logarithmic extinction coef¬ 
ficient for H/3/Ha. The observed reddening in magnitude 
E(B-V) was also calculated using the following relationship 
(Seaton 1979), 


c{Hp) = 0.4 XpE(B - V), (2) 

where the extinction parameter Xp = 3.615 (Fitzpatrick 
1999). The errors of the extinction c(H/l) and E(B-V) are 
calculated through standard error propagation of equations 
(1) and (2). We obtain c(H/3)=0.13±0.03, in agreement 
with the value quoted by WLB05 (c(H/3)=0.139) and E(B- 
V)=0.09±0.02. Our logarithm H/3 flux -11.56±0.03 is also in 
very good agreement with previous studies (-11.53, WLB05; 
-11.52, Bar87a) and the recently published catalogue of inte¬ 
grated Ha fluxes of Galactic PNe (Frew et al. 2013). Notice 
that the interstellar extinction c(H(3) of Vy 1-2, given by 
Frew et al. (2013) is significantly lower than the previous 
values. 

The nebular electron temperature (T e ) and density 
(N e ) are derived using the TEMDEN task in IRAF (Shaw 
& Dufour 1994). From the [N n] (6584+6548)/5755 and 
[O i] (6300+6363)/5577 diagnostic line ratios, we obtain 
T e = 10850+2100 K and 10400+5400 K, respectively, whilst 
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Figure 13. Observed spectra of Vy 1-2 taken with the 1.3—m telescope. Upper panel displays the spectrum with an exposure time of 
300s, covering the wavelength range of 4700-6800 A. Lower panels display two enlargements of the spectrum (left: 4700-5700 A, right: 
5700-6800 A) with an exposure time of 3600 sec, in order to highlight the weaker emission lines. 
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12 Evidence for a WR or WEL type binary nucleus in Vy 1-2 


Table 4. Electron temperature (in K) and electron density (in Table 5. Ionic and total ( e(X)=log(n(X)/n(H))+12) abundances 
cm -3 ) of Vy 1—2. derived from the observed fluxes of Vy 1—2. 



this work 

WLB05 

Bar78a 

McNabb 4 

T e ([0 III]) 

— 

10400 

9800 

— 

T e ([N II]) 

10850T2100 

- 

- 

- 

Te([0 II]) 

- 

16179 

- 

7943 

T e ([0 I]) 

10400T5400 

- 

- 

- 

N e ([S II]) 

4500A650 

1160 

5000 

- 

Ne([Cl III]) 

4330A1600 

- 

- 

- 

N e ([Ar iv]) 

4550±1300* 

3300 

- 

- 

Ne([Q II]) 

- 

4100 

- 

4466 


f T e [0 II] and N e [0 II] estimated by McNabb et al. are obtained 
using the O II optical recombination lines. 

| The contribution of He I 4712 A emission line to [Ar iv] 4711 A 
line, was calculated based on the work of Benjamin et al. (1999). 

N e are calculated as 4500±650 cnE 3 , 4330±1600 cm -3 
and 4550±1300 cm' 3 from the [S n] 6716/6731, [Cl hi] 
5517/5537 and [Ar iv] 4711/4740 diagnostic line ratios, re¬ 
spectively. No evidence of density variation with the ionisa¬ 
tion potential is found (see also Fig. 14). For direct compar¬ 
ison, Table 4 lists the N e and T e together with those derived 
in previous studies from collisionally excited and recombi¬ 
nation lines. 

Regarding the N e [Ar], it should be noticed that the 
[Ar iv] 4711 A line is blended with the He I 4712 A, 
[Ne iv] 4724 A and 4726 A lines. Therefore, to estimate 
the correct N e , the contribution of the latter to [Ar iv] 
4711 A was calculated and subtracted. The theoretical value 
of He i 4712 A/He i 4471 A is 0.146 for T e =10000 K, 
N e =10000 cm -3 and 0.105 for T e =10000 K, N e =1000 cm -3 
(Benjamin et al. 1999). The observed intensity of He I 
4471 A is 4.403 (relative to H^ = 100; WLB05). This value 
gives 1(4712) equal to 0.64 and 0.46 for N e =10000 cm -3 
and N e =1000 cm -3 , respectively. The observed intensity of 
He I 4712 A is 0.519 (WLB05) and lies between the the¬ 
oretical lower and upper limits, as it is expected, given 
that N e = 4330 cm -3 also lies between the equivalent lower 
(1000 cm -3 ) and upper (10000 cm -3 ) limits. Our best- 
fit photo-ionisation model predicts I(4712)=0.53 (see §4.2). 
The contribution of [Ne iv] was considered negligible since 
only a star with T e ff > 100 kK can significantly ionise this 
line. WLB05 found the [Ne iv]/H/3 equal to 0.216 or 5% of 
the [Ar iv] 4711 A line. 

The recently released software, namely PyNeb (Luridi- 
ana et al. 2011) was used to construct the T e -log(N e ) diag¬ 
nostic plot for Vy 1-2 from our data and those from WLB05 
(Fig. 14). This diagram illustrates the whole picture of the 
physical conditions in Vy 1-2 where a good solution for T e - 
log(Ne) is found; T e ~11000 K and log(N e ) ~ 3.6 cm -3 . 
Note that both high T e [0 n] (16200 K or 17110 K corrected 
for recombination excitation on the nebular lines) and low 
N e [S n] (1160 cm -3 ) derived by WLB05 do not agree with 
this solution. 

Most of the emission lines in our spectra are in very 
good agreement with those reported by WLB05, except the 
[S n] 6731 A emission line, which surprisingly differs by a 
factor of 1.5 (see Table 3) and reflects on the significant N e 
difference. Our estimation of N e ([S ll])=4500 cm -3 agrees 


ionic 

this work 

this work 

WLB05 

Bar78b 

abundances 

He+/H+ 

(ICF) 

0.087A0.008 

(model) 

0.083 

0.084 

0.076 

He++/H+ 

0.022T0.004 

0.024 

0.023 

0.019 

ICF(He) 

1.00 

1.00 

1.00 

1.0 

e(He) 

11.04T0.03 

11.03 

11.03 

10.98 

10 5 x 0+/H+ 

1.55i0.17“ 

2.35 

1.36 

2.51 

10 5 x 0 2 +/H+ 

37.74A2.78 

66.11 

41.55 

47.1 

ICF(O) 

1.17 

1.19 

1.18 

1.25 

e(O) 

8.66i0.04 

8.91 

8.70 

8.79 

10 5 x N+/H+ 

0.47T0.05 

0.62 

0.395 

0.435 

10 5 x N 2 +/H+ 

- 

- 

9.300 

- 

ICF(N) 

29.28 

41.45 

1.206 

24.8 

e(N) 

8.13T0.05 

8.43 

8.06 

8.03 

10 6 x S+/H+ 

0.24T0.04 

0.38 

0.144 

0.385 

10 6 x S 2 +/H+ 

2.71T0.37 

12.19 

2.467 

- 

ICF(S) 

2.17 

1.86 

2.334 

24.8 

<S) 

6.80T0.08 

7.37 

6.78 

6.98 

10 8 x Cl 2 +/H+ 

7.09A1.97 

10.07 

_ 

_ 

ICF(Cl) 

2.35 

2.03 

- 

- 

e(Cl) 

5.22T0.12 

5.31 

- 

- 

10 7 x Ar 2 +/H+ 

11.23il.35“ 

13.03 

9.370 

— 

10 7 x Ar 3 +/H+ 

6.26A1.43 

10.34 

6.902 

- 

10 7 x Ar 4 +/H+ 

0.75T0.26 

1.09 

1.005 

- 

ICF(Ar) 

1.04 

1.01 

1.035 

- 

e(Ar) 

6.28i0.06 

6.39 

6.25 

- 


“ The ionic abundances of these ions were calculated using the 
intensities from WLB05, considering an error of 10%, since no 
errors are given in WLB05. 

with the value reported by Baker (1978b) as well as the re¬ 
cent N e determinations from the O II optical recombination 
lines (ORLs) by Me Nabb et al. (2013). It is worth mention¬ 
ing the high difference of order 2 in T e derived from the O II 
optical recombination lines by WLB05 and McNabb et al. 
(2013). 

The ionic abundances of N, O, S, Cl and Ar, relative 
to H + , are derived using the ionic task in IRAF (Shaw & 
Dufour 19940. T e =10600 K and N e = 4500 cm 3 are used 
for both low and high-ionisation regions. 

The He + /H + and He ++ /H + ionic abundance ratios are 
calculated based on the work of Benjamin et al. (1999). For 
the total chemical abundances we used the ionisation cor¬ 
rection factor in order to correct for the unobserved states 
of ionisation (ICF; Kingsburgh & Barlow 1994). Since the 

3 Transition probabilities: N+, 0 + and O 24 ", Wiese et al. (1996); 
S 4 ", Verner et al. (1996); S 4 " 2 , Mendoza & Zeippen (1982); Cl 24 ', 
Mendoza (1983); Ar 44 ", Mendoza & Zeippen (1982a); Ar 44 ", Men¬ 
doza &c Zeippen (1982b). Collisional strengths: N 4 ", Lennon & 
Burke (1994); O 4 ", Pradhan (1976); 0 2+ , Lennon & Burke (1994); 
S+, Ramsbottom et al. (1996); S 4 " 2 , Galavis et al. (1995); Cl 24 ', 
Butler & Zeippen (1989); Ar 34 ", Zeippen et al. (1987); Ar 44 ", 
Galavis et al. (1995). 
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Figure 14. T e -log(N e ) diagnostic plot for Vy 1-2. 


aforementioned authors do not provide an ICF for the Cl, 
we used the equivalent equation from Liu et al. (2000). In 
Table 5, we list the observed ionic and total abundances de¬ 
rived in this work, together with those derived by WLB05 
and Bar78b. For a direct comparison the ionic abundances 
from our best-fit model are also included in Table 5. The 
ICFs computed from our model are listed in Table 5. One 
can see that N shows the largest discrepancy in ICFs be¬ 
tween the empirical method (Kingsburgh & Barlow 1994) 
and our photo-ionisation model. This difference in the ICFs 
is not such high to significantly alter the N/O abundance 
ratio and the classification of the nebula as non-Type I, like 
in NGC 1501 (Ercolano et al. 2004). 


3.3.2 Photo-ionisation modelling 

In order to perform a more assiduous study on the ionisa¬ 
tion structure of Vy 1 -2, we constructed a photo-ionisation 
model using Cloudy (version 13.02; Ferland et al. 2013). Sev¬ 
eral atomic, ionic and molecular databases are used within 
Cloudy (e.g. Dere at al. 1997; Badnell et al. 2003; Badnell 
2006, Landi et al. 2012; among others). The complete infor¬ 
mation can be found in Linkins et al. (2013) and Ferland et 
al. (2013). 

Cloudy code uses as input parameters, (a) the energy 
distribution of the ionising radiation, (b) the luminosity of 
the central star and (c) some additional constraints that 
are usually estimated from the observed spectrum, which 
should as well be reproduced by the model such as the T e 
and N e . Both parameters are required for determining the 
thermal balance of PNe. Because the configuration of the 
spectrograph does not include the lines of [O n] 3726 A and 
3729 A, [Ne hi] 3869 A and 3968 A, He n 4686 A, [Ar hi] 
7135 A, they were drawn from WLB05. 

The model of Vy 1 2 was developed using the simplest 
set of assumptions (i) a blackbody approximation for the 
central star and (ii) a spherically symmetric nebula with 
constant density. The initial abundances of He, N, O, S, Ar 
and Cl were taken from our empirical analysis whereas Ne 
and C abundances from WLB05, and all of them were kept 


as free parameters, except C. Dust grains were also included 
in the model since they play an important role in the thermal 
balance of PNe. However, there is no information regarding 
the grain composition (no ISO spectra are available) and a 
mixture of Polycyclic Aromatic Hydrocarbons (PAHs; opti¬ 
cal constants from Desert et al. 1990, Schutte et al. 1993, Li 
& Draine 2001; for more detail see Ferland et al. 2013) and 
astronomical graphite and silicate (Martin & Rouleau 1991; 
Laor & Draine 1993) with a constant dust-to-gas ratio, was 
adopted. The amount of PAH, graphite and silicate, grains 
were also kept as free parameters. Furthermore, the filling 
factor was kept as a free parameter in order to reproduce 
the observed H/3 flux. 

Given that C, O and N are the most important coolant 
elements in PNe, and they play a significant role in the ther¬ 
mal balance, small changes in their abundances can signif¬ 
icantly affect the intensities of the emission lines of all ele¬ 
ments. Hence, in order to get the best fit between the mod¬ 
elled and observed spectrum, we first varied the N and O 
abundances until the intensities of [N n] 6584 A and [O ill] 
5007 A emission lines were fitted. Then, regarding the other 
elements, we varied the abundances of the rest of elements 
(He, S, Ar and Cl) up to the point that the intensities of 
He ii 5876 A, [S n] 6716, 6731 A, [Ar iv] 4712, 4740 A and 
[Cl ill] 5517, 5538 A lines, were also fitted. The modelled 
spectrum of Vy 1-2 is presented in Table 3 together with 
the observed spectrum. The mean difference between the 
observed and modelled intensities is ^10% for the bright 
lines (S/N>5) and ~20-30% for the fainter lines (S/N< 5). 
The [K iv] 6101 A and [Si Ii] 5040 A lines, as well as the 
stellar C II 6461 A, C iv 5801 and 5812 A and the O ill 
5592 A lines, were not used to constrain the model. 

As for the carbon, we assumed a constant abundance of 
log(C/H) + 12=8.08, which was derived from the collisionally 
excited lines but it was not possible to get a good match of 
the CEL and ORL at the same time (very weak ORL), while 
the ionisation balance of C was far from the observed value. 
Assuming a higher C abundance (of 9.0), the model better 
predicts the ORL but overestimates the CEL. In addition 
to the blackbody stellar model, a H-deficient stellar atmo¬ 
sphere model for the ionizing source was also used. A good 
match between the modelled and observed spectrum could 
be obtained only for the case of high C abundance (e(C)>9). 

In conclusion, whatever stellar atmosphere model was 
chosen for Vy 1-2, we did not fully reproduce the observed 
properties (e.g. the ionisation balance of the elements) to a 
satisfactory degree. In particular, the ionisation balance of 
all the elements is overestimated by the low C abundance 
model, whilst the higher C abundance model seems to pro¬ 
duce an ionisation balances that matches better the observed 
ones. Also, our simple spherically symmetric shell geometry 
seems to be far from the real structure of the nebula, which 
as discussed in § 3.1 seems to have a bipolar structure or 
at least a significant deviation from a spherically symmetric 
morphology. 

4 DISCUSSION 

4.1 Morpho—kinematic analysis 

Vy 1 -2 shows a high-speed component expanding with ve¬ 
locities up to 90-100 km s -1 , which is likely associated with 
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Table 6. Total abundances, e(X)=log(n(X)/n(H))+12, of Vy 1-2. 


6(X) 

this work 
(ICF) 

this work 
(model) 

WLB05 

Bar78b 

Stan06 

Solar a 

He 

11.04 

11.03 

11.03 

10.98 

11.23 

10.93 

C 

- 

8.08 6 

8.08 6 

- 

- 

8.43 

O 

8.66 

8.91 

8.70 

8.79 

8.46 

8.69 

N 

8.13 

8.42 

8.06 

8.03 

7.90 

7.83 

S 

6.80 

7.37 

6.78 

6.98 

- 

7.12 

Ne 

- 

8.11 

8.00 

8.19 

7.81 

7.93 

Ar 

6.28 

6.39 

6.25 

- 

5.84 

6.40 

Cl 

5.22 

5.31 

- 

- 

- 

5.50 

log(N/0) 

-0.53 

-0.49 

-0.64 

-0.76 

-0.56 

-0.86 

l°g(S/0) 

-1.86 

-1.54 

-1.92 

-1.81 

- 

-1.57 

log(Ne/0) 

- 

-0.80 

-0.70 

-0.60 

-0.65 

-0.76 

log(Ar/0) 

-2.38 

-2.52 

-2.45 

- 

-2.62 

-2.29 

l°g(Cl/0) 

-3.44 

-3.60 

- 

- 

- 

-3.19 


a Asplund et al. 2009 

b C abundance was derived from collisionally excited lines. A higher C abundance of e(C)=8.98 was also derived from the optical 
recombination lines (WLB05). Vy 1-2 is a nebula with moderate high abundance discrepancy factors. 


the material in the faint bipolar lobes, with the eastern and 
western being blue- and red-shifted, respectively. Consider¬ 
ing that these lobes are inclined by 10° with respect to the 
line of sight and their projected distance on the sky is 3.6 " 
(measured from the eastern to the western lobe), we get a 
rough estimation of the nebula’s size along the polar axis 
of 20.4" (r=10.2"). Given that the bipolar lobes have a de- 
projected radial expansion velocity of 90-100 km s -1 , their 
kinematic age is found to be 5200±750 years, for D=9.7 kpc. 

A pair of knots is marginally detected in the HST [O ill] 
image, located at 11.5" (knot S) and 14.5" (knot N) from 
the central star. The projected expansion velocity of knot N 
differs from the I4 ys of the nebula by 5 km s _1 . This converts 
to a de-projected radial velocity in a direction away from the 
central star of V=5/cos($) km s -1 , where 9 is the inclination 
angle of the knot to the line of sight. If we consider that the 
pair of knots is the result from the same ejection that formed 
the bipolar nebula (same age), its de-projected expansion 
velocity should be up to ~130 km s -1 . Taking into account 
the projected velocity of knot N (5 km s -1 ), as measured 
from the PV diagrams (Fig. 8), the inclination angle of the 
knots is found 89° with the respect to the line of sight. This 
suggests a possible precession of the putative central star 
system and multiple ejection events leading to the formation 
of a bipolar nebula. 

A single Gaussian was fitted to the emission line profiles 
of Ha, He II, [N n] and [O ill] and gave expansion veloci¬ 
ties for the unresolved ring-like structure of 19±3 km s -1 , 
(23-25)±4 km s -1 , (8-9)±2 km s -1 and 18±3 km s -1 , re¬ 
spectively. For Vex P =19 km s _1 and r=1.4", we estimate the 
kinematic age of the ring-like structure ~3500±500 years. 
This means that the ring-like structure was formed more 
recently than the bipolar lobes, which is not feasible. If, al¬ 
though, the inclination angle of the bipolar lobes is slightly 
higher of 15°, we get the same dynamical age for both com¬ 
ponent, ring-like and bipolar lobes, of 3500 years. 

The velocity trend, on which the low-ionisation gas 
(e.g. [N il]) expand slower than the high-ionisation gas-e.g. 


[O ill], He II, is not consistent with a typical homologous 
expansion law, where their expansion velocity increases lin¬ 
early with the distance from the central star. Despite the 
fact that V[o hi] and VHe ii could be roughly considered 
the same within the uncertainties, the much lower Vr N n i 
provides strong evidence for a non-homologous expansion. 
A small group of PNe (e.g. BD +30°3639, Akras & Stef¬ 
fen 2012; see also Medina et al. 2006) has been found to 
show this velocity trend. The central stars of those PNe are 
either [WR] or WEL-type stars (Medina et al. 2006). The 
acceleration of the inner zones of the nebula might be the 
result of strong, fast winds from a [WR] or WEL nucleus 
(Gesicki & Zijlstra 2003). The detection of the C II 6461 A, 
C IV 5801 and 5812 A and O ill 5592 A stellar lines in 
Vy 1-2 supports this scenario (see §4.3.2). Another possible 
scenario responsible for the acceleration of the inner nebu¬ 
lar regions is a nova-type eruption of a close binary system. 
Akras & Steffen (2012) propose this mechanism in order to 
explain the morpho-kinematic characteristics of the young 
BD +30° 3639 nebula. 


4.2 Chemical analysis 

The plasma parameters of Vy 1-2 such as T e , N e and chem¬ 
ical abundances were derived by means of the empirical 
methods and photo-ionisation modelling. Our best-fitting 
model provides a good approximation of the observed spec¬ 
trum, except for the [O hi] 4959, 5007 A doublet, and adds 
some constraints to the nebular and stellar properties of 
Vy 1-2. The predicted H/3 flux is in a good agreement with 
the observed flux (see Table 3) with a filling factor equal to 
0.44. The mean difference between the observed and mod¬ 
elled lines varies from 10% for the bright lines (S/N> 5; 
the discrepancy in the [O ill] 4959, 5007 A doublet is about 
14%) to ~20-30% for the fainter lines, S/N<5. Moreover 
the intensities of [O i] and [N i] neutral lines are underesti¬ 
mated. There are two possible scenarios for this: i) there is a 
large amount of neutral material around the nebula and/or 
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Figure 15. Hertzprung-Russell diagrams for hydrogen-burning (left panels) and helium—burning (right panels) models. Upper panels 
display the evolutionary tracks from Vassiliadis &; Wood (1994) produced by (Mi n j t i a i, Mfi na i) = (IMq, 0.569Mq), (1.5Mq, 0.597Mq) 
and (2Mq, 0.633Mq) for the hydrogen—burning model (left panel) and (IMq, 0.567Mq), (1.5Mq, 0.600Mq) and (2Mq, 0.534Mq) for 
the helium-burning model (right panel). Lower panels display the evolutionary tracks from Blocker (1995) produced by (M; n ; t ; a i, Mf; na i) 
= (3Mq, 0.605Mq), (3Mq, 0.625Mq) and (4Mq, 0.696Mq) for the hydrogen-burning model (left panel) and (IMq, 0.524Mq), (3Mq, 
0.625Mq) and (5Mq, 0.836Mq) for the helium—burning model (right panel). The position of the central star of Vy 1-2 for a distance of 
9.7 kpc is marked with a red star symbol. 


ii) shock excitation of neutral lines close to the ionisation 
front may be significant. We can actually rule out the first 
scenario since the model predicts a density-bounded neb¬ 
ula whereas high expansion velocities up to 90-100 km s -1 
could favour the second scenario. 

The total chemical abundances of Vy 1-2, for a black- 
body stellar flux distribution, are listed in Table 6, which 
also includes those derived from previous studies, as well 
as the solar abundances, for direct comparison. The low He 
and N/O abundances imply a non Type-I nebula based on 
both Peimbert & Torres- Peimbert’s (1983) and Kingsburgh 
& Barlow (1994) classifications. The O abundance derived 
from the empirical ICF is found to be almost solar and close 
to the average value of non-Type I PNe (Bohigas 2008). Our 
best-fitting model predicts a somewhat higher value by 0.25 
dex meaning that the star has not undergone the ON-cycle 
burning during the second dredge-up phase. 

Given that Vy 1-2 is an O-rich nebula - log(C/0)=- 
0.62, WLB05; or log(C/O)=-0.83; our model - and the C 
abundance (e(C)=8.08) is lower than half of the solar (As- 


plund et al. 2009) as well as the average non-type I PNe val¬ 
ues, the N enrichment must have occurred by the conversion 
of C to N, via the hot-bottom burning (HBB; CN-cycle). It 
is known from theoretical models of stellar nucleosynthesis 
that only the progenitor stars with masses ^3 M@ undergo 
the 2 nd dredge-up and HBB events (Karakas et al. 2009). 
This result is inconsistent with the low N/O abundance ra¬ 
tio (-0.53 (ICF), -0.48 (model)) and a low-mass progenitor 
star (^1 Mq). 

Unlike the single star evolutionary models, a close bi¬ 
nary system can give such low N/O and C/O abundance ra¬ 
tios via the mass-transfer exchange (Iben & Tutukov 1993, 
De Marco 2009). Hf 2-2, NGC 2346 and Hb lfl are PNe 
that have confirmed close binary systems and low N/O and 
C/O values (see De Marco 2009). In a common envelope 
evolution, the AGB-phase may be significantly shorter than 

4 Despite that Hsia et al. (2006) presented the detection of pe¬ 
riodic variabilities, the binarity of the central star of Hb 12 still 
remains under debate (De Marco et al. 2008) 
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in a single star case resulting in the formation of an O-rich 
envelope (Izzard et al. 2006). 

A comparison of the chemical abundance of Vy 1-2 with 
those derived from PNe with different types of nucleus (nor¬ 
mal, WEL and [WR]) shows that the He abundance and 
N/O ratio better match with those from WEL-type nuclei 
while the O abundance with those from WEL- and [WR] 
types (Corny et al. 2009). As for the other elements, Ne 
matches those from WEL-type nuclei, Ar those from WEL- 
and normal-type nuclei, while Cl abundance is found to 
be significantly lower than those values derived by Corny 
et al. (2009) for all types of nuclei, but closer to the so¬ 
lar value. In particular, the average abundance of Cl varies 
from 6.12 to 6.48, while we find e(Cl)=5.22 (ICF method) 
and 5.31 (photo-ionisation model). S abundance seems to 
match all types of nuclei due to the large difference between 
the ICF method (e(S)=6.80) and the photo-ionisation model 
(e(S)=7.37). This may be a consequence of the poor fit of 
the [S ill] 6312 A emission line. The abundance uncertain¬ 
ties are ~10% for the He and ~20-30% for the rest of the 
elements. 

It is worth mentioning here that Vy 1 -2 has also been 
considered as a candidate halo PN due to its large distance 
from the Galactic plane (~500) and its high systemic veloc¬ 
ity (-85+4 km s _1 ) (Quireza et al. 2007). However, a com¬ 
parison of its chemical abundances with those from Galactic 
disk and halo PNe (Otsuka et al. 2009) does not support this 
classification. 

The best-fitting model of Vy 1-2, presented here, fails 
to reproduce the [O ill] 5007 A and 4959 A emission lines. 
The difference between the observed and modelled values is 
about 14%, similar to another PN with a WR-type central 
star as well as large abundance discrepancy factors (ADFs), 
NGC 1501 (Ercolano et al. 2004). We have to point out, 
though, that the emission line ratio [O ill] 5007/4363, which 
is an indicator of T e of the nebula, is predicted with an 
accuracy better than 5%. The main causes of the model’s 
inability to reproduce the observed ionization balance and 
predict the bright [O ill] lines may be: i) the structure of the 
nebula, which is far from the spherically symmetric assump¬ 
tion, ii) the stellar continuum, where neither the blackbody 
nor an H-deficient stellar atmosphere model were able to 
reproduce all the observed data and/or iii) the dust compo¬ 
sition, which seems to play a crucial role in this objects, but 
no information is available for Vy 1-2. 

The large ADFs found in Vy 1-2 for 0 ++ , N ++ , C 3+ , 
C++ and Ne++ are of 6.17, 11.8, 14.0, 9.27 and 13.8, respec¬ 
tively (see WLB05) and could indicate the possible presence 
of cold H-deficient clumps in the nebula, which may also 
hinder the results of our model This possibility is, however, 
beyond the scope of this paper. 


4.3 Central star 

4-3.1 Evolutionary status 

Our simple blackbody stellar model predicts a luminosity 
log(L/L0)=3.4O and an effective temperature T e fi=98 kK, 
for D=9.7 kpc. For a distance of 7.8 kpc, we obtain 
log(L/LQ)=3.21 in good agreement with those values re¬ 
ported by Stangellini et al. (2002) and Phillips (2003). Fig¬ 
ure 15 portrays the positions of Vy 1-2’s central star in 


the HR diagram, for a distance of 9.7 kpc. The theoretical 
evolutionary tracks plotted are the hydrogen-burning (left 
panel) and helium-burning (right panel) models from Vassil- 
iadis & Wood (1994; upper panels) and (Blocker 1995; lower 
panels). In Vassiliadis and Wood’s models, the central star 
corresponds to a progenitor mass of 1 M@ and a core mass 
of 0.569 M 0 (hydrogen-burning ) and 0.567 M Q (helium¬ 
burning). The timecales from these evolutionary tracks are 
significantly higher than the kinematic age of the nebula 
with 26000 years and 41000 years for the hydrogen- and 
helium-burning respectively. 

In Blocker’s models, the evolutionary track for a very 
late thermal pulse (VLTP) (lower, right panel) implies a 
progenitor star of 3 M@ and a core mass of 0.625 M@. Even 
though the timescale of the former evolutionary track is still 
higher than the kinematic age of the nebula, by a factor 
of ^2, it is closer to the observed value. The scenario of 
a close companion and mass-transfer exchange would be 
able to explain the rapid evolution of the star. Unlike the 
1 M@ low mass star, a 3 M@ progenitor star can go through 
the HBB effect during the 2 nd dredge-up phase resulting in 
enhancement of N and depletion of C. The detection of the 
C II 6461 A, the doublet C ivA A 5801, 5812 and the O ill 
5592 A stellar lines emitted from a hydrogen-deficient star 
supports the VLTP scenario. 

Our new distance value indicates that Vy 1 -2 is a dis¬ 
tant nebula (9.7 kpc) and provides a much better agreement 
between the kinematic age of the nebula and the evolution¬ 
ary age of the central star. 

4-3.2 [WR] or WEL type nucleus? 

Apart from the detection of the common nebular emission 
lines, weak stellar emission lines such as C II 6461 A, C IV 
5801 and 5812 A and O ill 5592 A were also detected for 
the first time in this nebula in our very deep spectra (Table 
7). Their line profiles are shown in Fig. 16. There are two 
groups of stars with emission lines of highly ionised He, C, O 
and N: [WR] (Crowther 2008) and WELS stars (Tylenda et 
al. 1993). Although, both classes exhibit the same emission 
lines (e.g C IV doublet 5801 and 5812 A and [C ill] 5696 A), 
the former show much broader lines profiles and stronger 
intensities than the latter due to the fast and strong stellar 
winds. 

Based on the classification scheme of Acker & Neiner 
(2003) and the C n/C IV and O ill/C IV line ratios, the cen¬ 
tral star of Vy 1-2 is classified as a late-type star ([WRL]). 
On the other hand, the absence of the [C ill] 5696 A line (see 
WLB05) in conjunction with the detection of the C IV dou¬ 
blet, N ill and C III-IV 4650 A lines, also suggest a WEL- 
type nucleus. Therefore, new high -resolution observations 
are required in order to confirm the spectral type of this 
star. 


4.4 Infrared properties and SED 

Vy 1-2 is also a strong infrared emitter indicating the pres¬ 
ence of significant amount of dust. Thus far, no one had 
studied the infrared emission and dust properties of this 
compact nebula. For this study, we used infrared data from 
2MASS, IRAS (IRAS 17524+2800), AKARI and WISE data 
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Figure 16. Part of the spectrum of Vy 1-2, in which the faint 
stellar lines C iv A A 5801, 5812 A, C n 6461 A and O m 5592 A 
are seen (in units of 10 -14 erg s ’ 1 cm -2 ). 


archives. The data cover a wavelength range from the near- 
to the far-infrared. 

The spectral energy distribution (SED) of Vy 1-2 is pre¬ 
sented in Fig. 17 and displays a double-peaked profile. One 
peak lies at the optical wavelength due to the hot central 
star (100,000 K) of the nebula and the second one at 30- 
40 /im implying the presence of a cold dusty circumstellar 
shell with a temperature of ^100-150 K. Vy 1-2 may also 
exhibit a small K-band excess, associated either with a hot 


Table 7. Stellar emission lines of Vy 1-2 


Line (A) 

F(A)“ 

S/N 

O III 5592 

0.08 

3.0 

C iv 5801 

0.12 

3.0 

C iv 5812 

0.11 

3.0 

C ii 6461 

0.09 

3.0 


“ Relative to the nebular H/3=100 


dust component or the emission from a cool (~2000 K) un¬ 
resolved companion. 

The all-sky Wide-Field Infrared Survey Explorer 
(WISE) observed Vy 1-2 in four photometric bands at 
3.35 /im, 4.6 yum, 11.6 /im and 22.1 yim (W1 through W4, re¬ 
spectively; Fig. 18). Because of the low angular resolution in 
these bands (6.1", 6.4", 6.5"and 12", respectively; Wright et 
al. 2010), Vy 1-2 remains unresolved. It is worth mentioning 
here that W1 and W2-band images unveil the presence of 
a seemingly collimated outflow in the same direction as the 
NW knot (PA=310°). An arrow in Figure 18 indicates the 
position of the collimated structure. Despite the W3-band 
having a similar angular resolution to those of W1 and W2 
bands, the collimated structure can not be discerned. 

By scrutinizing the DSS images for Vy 1-2, we discover 
the presence of a very faint starlike object in the same di¬ 
rection as the outflow and in the same position as the NW 
knot (~14" away from the central star). The detection of 
this held star in the DSS IR band, as well as in the W1 and 
W2-bands but not in the W3-band, implies a very cold star 
with a temperature of ~2000 K. Note that the seemingly 
halo observed in the W3 and W4 band images (panels (c) 
and (d) in Fig. 18) is not a real structure since it is also 
apparent around the bright held stars. 

WISE bands, encompass the dust continuum and sev¬ 
eral features attributed to PAH features at 3.3 /im, 8.6 /im, 
11.3 /im, 12.7 /im and 16.4 /im, silicate dust features at 
9.7 /im and 18 /im as well as high excitation nebular lines 
(e.g. [S vi], [Ne iv] and [Ne v]). However, ISO spectra are not 
available, prohibiting any further investigation on the dust 
composition. Nevertheless, the detection of single ionised re¬ 
combination C lines (C n 6578 A and 6581 A) may suggest 
the presence of C-rich dust. The large excess in the W3 and 
W4 bands may also be associated with the presence of sili¬ 
cate dust features. The SwSt 1 nebula exhibits these C lines 
as well as the silicate dust feature at 10 /im (De Marco et al. 
2001). A double-dust chemistry (O- and C-rich) is observed 
in PNe with [WR] or WEL type central stars, probably due 
to the VLTP that turn the central stars from O-rich to C- 
rich. 

In order to constrain the dust emission, we used the 
aforementioned broad-band infrared huxes. The best-fitting 
SED to the infrared data was possible only assuming a 
C-rich dust model. Figure 17 displays the modelled SED 
of Vy 1-2. The PAHs are ~5.7 times more abundant (by 
mass) than the silicates and the dust-to-gas ratio equal to 
1.05xl0~ 2 . The total mass of the nebula and the dust mass 
are 0.08 M0 and 8.4 10 -4 M©, respectively. The high uncer¬ 
tainty of Vy 1-2’s distance implies as well high uncertainties 
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Figure 17. The observed and modelled SED of Vy 1-2 in the 
wavelength range from 0.43 to 140 /im. The data are from GSC2.2, 
GSC2.3, USNOB1, USNOB2, 2MASS, AKARI and WISE sur¬ 
veys. 

in these masses. This result is consistent with the scenario of 
a hydrogen deficient star and dual-dust chemistry. Interest¬ 
ingly, the ionized gas indicates an O-rich nebula (C/0<1), 
whilst the dust composition is C-rich. One possible explana¬ 
tion, as recently proposed by Guzman-Ramirez et al. (2014), 
is the formation of PAHs from the photo-dissociation of CO 
in a torus, such as the bright ring-like structure in Vy 1-2, 
whereas the dust in the main nebula remains O-rich. More¬ 
over, a recent change of the stellar surface chemistry due 
to the second dredge-up of C, as has been proposed to ex¬ 
plain the dual-dust chemistry in the young PN BD +3036 
(Guzman-Ramirez et al. 2015), can also provide a plausi¬ 
ble explanation for Vy 1-2 as well. We have to stress out 
though there is no available infrared spectroscopic data for 
Vy 1-2 and our analysis is based on the results from our 
best-fitting photo-ionisation model using the broad-band 
infrared fluxes. 

In summary, a careful look at the modelled SED reveals 
a lower flux in the infrared H-band, indicating that this mea¬ 
surement may be underestimated and the seemingly K-band 
excess is not really an excess associated with the presence 
of a cold companion as previously mentioned. 


5 CONCLUSIONS 

The main conclusions of this work are summarized as follow: 

1) Vy 1 -2 presents a bipolar structure, seen al¬ 
most pole-on. The bipolar lobes expand with velocities 
up to 90-100 km s -1 , whereas the bright, inner ring¬ 
like structure shows lower expansion velocities of Vh 0 
= 19 km s -1 , VHe ii =23-25 km s -1 , V[ N n]=8-9 km s -1 
and V[o m]=18 km s -1 . These velocities indicates a possi¬ 
ble acceleration process of the inner nebular regions, which 
may be associated with the strong/fast stellar wind from a 
[WR] or WEL nucleus. The kinematic age of the nebula is 
~3500 years for a distance of 9.7 kpc. 

2) A pair of knots is also apparent along P.A.=305°, 
located at 11.5" and 14. 5" and moving perpendicular to the 
line of sight. Assuming that they were formed approximately 
at the same time as the nebula, we estimated their velocities 
up to a value of 130+15 km s -1 . 

3) Vy 1-2 is apparently an O-rich, non-type I 



Figure 18. WISE images of Vy 1-2 in four bands: a) 3.35 (im, 
b) 4.6 /im, c) 11.6 /im and d) 22.1 fin i. The arrows indicate the 
position of the field star. The field of view is 180” X ISO”. North 
is up and East is to the left. 


nebula (log(C/O)=-0.58 (ICF) and -0.83 (model), - 
0.53^1og(N/0)^-0.49). Compared to solar abundances, O 
was found to be solar, C subsolar and it is enriched in N. 
The HBB CN-cycle must have occurred in order to explain 
these abundances, implying a progenitor star with initial 
mass ^3 Mq. 

4) None of the stellar ionizing continua used in our 
photo-ionisation model (blackbody or H-deficient) were 
able to reproduce the observed ionisation balance of the el¬ 
ements. The spherical symmetry asssumed for the nebula 
may be far from the real structure of Vy 1-2, and a more 
sophisticated 3-D bi-abundance model is required for the 
proper study of its ionisation structure. 

5) The low C/O and N/O abundance ratios of Vy 1- 
2 can not be attained by single star evolutionary models, 
and a post-common envelope binary system seems a more 
plausible scenario. The N enhancement can be explained by 
the binary interaction via mass-transfer exchange. 

6) A very late thermal pulse may also have occurred 
resulting in the formation of a H-deficient star. The detec¬ 
tion of the C ii 6461 A, the C ivA A 5801, 5812 doublet and 
the O ill 5592 A stellar lines supports this scenario. The 
timescale of VLTP evolutionary track for a 3 Mq progenitor 
star is consistent with the kinematic age of the nebula. 

7) Vy 1-2 is found to be very bright at the mid-IR 
indicating large amount of dust. Its SED shows a double- 
peaked profile with one peak at 30-40 fim which implies a 
circumstellar envelope with a temperature of ~100 K. The 
C/0<1 abundance ratio of the ionised gas indicates a O- 
rich nebula, whilst our best-fitting model predicts a C-rich 
dust, with PAHs being almost 5 times more abundant (by 
mass) than silicates. A possible interpretation of this result 
may be the formation of PAHs in a torus due to the photo¬ 
dissociation of CO and the formation of silicates in the main 
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nebula or a recent change of stellar chemical composition due 
to the second dredge-up of C. 
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